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Abstract. Using Hankel operators and shift-invariant subspaces on Hilbert space, this pa- 
per develops the theory of the operators associated with soft and hard edges of eigenvalue 
distributions of random matrices. Tracy and Widom introduced a projection operator W 
to describe the soft edge of the spectrum of the Gaussian unitary ensemble. The subspace 
WL 2 is simply invariant under the translation semigroup e ttD (t > 0) and invariant under 
the Schrodinger semigroup e lt ( D (t > 0); these properties characterize WL 2 via Beurl- 
ing's theorem. The Jacobi ensemble of random matrices has positive eigenvalues which 
tend to accumulate near to the hard edge at zero. This paper identifies a pair of unitary 
groups that satisfy the von Neumann- Weyl anti-commutation relations and leave invari- 
ant certain subspaces of L 2 (0, oo) which are invariant for operators with Jacobi kernels. 
Such Tracy-Widom operators are reproducing kernels for weighted Hardy spaces, known 
as Sonine spaces. Periodic solutions of Hill's equation give a new family of Tracy-Widom 
type operators. 
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1 Introduction 

This paper concerns the spectral theory and invariant subspaces of operators that arise 
in random matrix theory, particularly the soft and hard edges that occur on the limiting 
eigenvalue distributions of the Gaussian and Jacobi unitary ensembles. Tracy and Widom 
[28, 29, 30] introduced various operators to describe the soft edge of the spectrum of the 
Gaussian unitary ensemble; that is, the eigenvalues near to the supremum of the support 
of the equilibrium distribution. Burnol proposed that the theory of random matrices 
should be expressed in terms of Sonine spaces [7, p 692]. Here we develop this theory in 
a systematic manner to show that the Tracy-Widom calculations are instances of more 
general results on Hankel operators, and introduce new settings where the theory applies. 
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In section 2 we consider operators on L (R) with kernels 

w( X , v) =^m^m (L1) 

where col [A, B] satisfies a first-order linear differential equation, and give sufficient con- 
ditions for W to be the square of a Hankel integral operator. Further, we show that the 
determinants det(J — zW) are related to the solutions of Marchenko integral equations. As 
we show in section 3, kernels such as W arise as reproducing kernels for weighted Hardy 
spaces on the upper half-plane C+ = {z : Qz > 0} as in [2, 6]. 

Let J be the 'flip' map Jf(x) = f(—x), and M u the multiplication operator / i— > uf. 
The classical Hardy space H 2 consists of the holomorphic functions F on C+ such that 
sup y>0 I^oo \F( X + w)\ 2 dx < oo, and we identify such a function with its L 2 boundary 
values. The Fourier transform is •?-"/(£) = e~ lx ^f(x)dx/\ /f 27r. Given u G L°°, the 
bounded linear operator \f 7 kKT*M u T* is the Hankel operator T u on L 2 (R + ) that has 
distributional kernel T*u{x + y) as in [23]. 

We recall how examples of such operators appear in the theory of the Gaussian unitary 
ensemble. Let xj t k and yj t k (1 < j < k < n) be a family of mutually independent N(0, 1/n) 
random variables. We let X n be the n x n Hermitian matrix that has entries [^n]jfc = 
(xj,k + iyj,k)/V2 for j < k, [X n ]jj = xjj for 1 < j < n and [X n ] kj = (x jjk - iyj,k)/V^ for 

(2) 

j < k; the space of all such matrices with the probability measure a n forms the Gaussian 
unitary ensemble. 

Bulk of the spectrum. The eigenvalues of X n are real and may be ordered as Ai < 



< A n , so their positions are specified by the empirical distribution fi n = (1/n) $^?=i $, 



A, 



As n — > oo, the empirical distributions converge weakly to the equilibrium distribution, 
namely the Wigner semicircle law 

p(dx) = — l { - 2 ,2]{x)y/A - x 2 dx, (1.2) 

(2) 

for almost all sequences (X n ) of matrices under <g>'£L 1 crn ' . The bulk of the spectrum 
consists of those eigenvalues in [—2, 2]. See [22, p 93]. 
Let B t be the operator on L 2 (TV) that has kernel 

smtn(x-y) 

B t (x,y) = — , 1.3 

n(x - y) 

let Is be the indicator function of a set S, and let P{ a ,(3) be the orthogonal projection 
on L 2 (R) given by P( a ,(3)f(x) = I( a ,^)(x)/(x); for brevity we write P+ = P(o,oo) and 
P_ =P ( _ OO)0) - 
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Let E an (k; a, (3) be the probability with respect to that (a, j3) includes exactly k 
eigenvalues. Mehta and Gaudin [22] showed that 



k 



E an (k;a,P) - LA. (^) z=1 det[;-^ M BiP (o ,4 (1.4) 

This determinant can alternatively be expressed in terms of the operator 
^ a : L 2 [-a,a] — ► L 2 that has kernel ^ a (x,y) = e lxy l[_ a:a ](y)/V27v and satisfies ^ a ^l = 

B a /ir- 

Hard edges. Let Y n random n x n matrices with independent N(0,l/n) entries, 
and let < Ai < . . . < A n be the eigenvalues of the positive operator Y*Y n . Then 
v n = ^ Y^j=i is the empirical eigenvalue distribution, and v n converges weakly almost 
surely to the Marchenko-Pastur distribution; so that 



f(x)u n (dx)^ y o m]/—^ (i.5) 

almost surely as n — > oo for all continuous and bounded real functions /. Thus the Xj 
tend to accumulate near to their minimum possible value of zero, where the density of the 
limiting distribution is unbounded; this is the hard edge effect. 

Hard edges also arise from random matrices of the Jacobi and Laguerre ensembles. 
The Jacobi ensemble of order N with parameters u, 7 > —1/2 at inverse temperature (3 > 
is the joint distribution function 

1 N 

°N } j(dx) = -z- + ^') /37 ( 1 - x if V 1 1 (** - *if dx i • • • dx N (1.6) 

JV j = l l<j<k<N 

where — 1<^i<...<xtv<1 are the eigenvalues. 

Forrester [12] showed that the integral operator F a,b on £ 2 ((0, l),dx) with kernel 



F - (x, y) = I (a>6) (x) * 2(x-y) — (y) (1-7) 

determines the limiting distribution of scaled eigenvalues Xj/(4n) from the Laguerre en- 
semble near to the hard edge, and conjectured that a similar result holds for the Jacobi 
ensemble. Using the orthogonal polynomial technique, Forrester and Rains [13] have veri- 
fied the cases of (3 = 1, 2 and 4 following earlier work by Borodin [5] and Duenez. 

We introduce the scaled eigenvalues £j by Xj = cos^/v^iV, to ensure that the mean 
spacing of the £j is of order O(l) near to the hard edge at Xj ~ 1. One can show that 

<T$ tJ [(a, b) contains no -> det(J - F a > b ) (AT -> 00). (1.8) 
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For subsequent analysis we change variables by writing x = e~ 2? and y = e~ 2v so that 
£, i] E (0, oo) for x, y G (0, 1). Let Gt be the unitary integral operator on L 2 (R) that has 
kernel e - ^ -7 ' J w (e~ i ~^~ v ); let Qi = GiP + Gi (£ G R), which gives a strongly continuous 
family of orthogonal projections. For compact operators S and T on Hilbert space, the 
spectrum of ST equals the spectrum of TS. The integral operator <E>^ = P + GeP+ on 
L 2 (0, oo) is Hilbert-Schmidt, and when < a < 1 and a = —(1/2) loga satisfies 

det(J - zF°' a ) = det(J - z$ 2 (a) ). (1.9) 

In section 4 we interpret these operators on the Sonine spaces u v H 2 where u v (x) = 
2 ix T((l + v + ix)/2)/r((l + v - ix)/2). 

Soft edge of the spectrum. We recall some results of Tracy and Widom [12, 28] con- 
cerning the largest few eigenvalues. The Airy function Ai(x), as defined by the oscillatory 
integral 

Ai(z) = — / e ^ zt+t3 ^ dt, (1.10) 

27T J-oo 

satisfies the Airy differential equation [27, page 18] y" — xy = 0. Let be the integral 

operator on L 2 (R) defined by the Airy kernel 

W 1/3 (z,v)= M(xWiy) - Ai ' iX)Aiiy) - (LID 
' x — y 

We scale the eigenvalues of the Gaussian ensemble by introducing 

= n 2/3 (^L,-2), (1 . 12) 



(2) 

and let E an (k; a, (3) be the probability with respect to oVi that (a, (3) contains exactly 
k of the £j (j = l,...,n); see [22, page 116, A7]. Aubrun [3] proved that the operator 
Wy'P = P( a ,p)Wi/3P( aj p) on L 2 (R_|_) is of trace class for < a < (3 < oo, and 

I ~\\k , jk 

E an (k-,Z-,a,P)^^(^) z=i det(l-zWffl (n - oo). (1.13) 

The compression of to L 2 (a, oo) may be identified, under the change of variables 

s^a + s, with r 2 ^ where the Hankel integral operator r( a ) on L 2 [0, oo) satisfies 

P oo 

T (a) f(s)= Ai(a + s + t)f(t)dt (/GL 2 (0,oo)). (1.14) 
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The spectrum of P( a ,/3)^fo-)P{a,i3) equals the spectrum of T^P^ a ^T^, and hence 

det(J - zP {ctj0o) W 1/3 P {aj00) ) = det(J - zT 2 a) ). (1.15) 

Edge distributions and KdV. For < t < 1 let w(x;t) be the unique solution to the 
Painleve II equation w" = 2w 3 + xw that satisfies w(x;t) x — ^/iAi(x) as x — > oo. By the 
theory of inverse scattering for the concentric Korteweg-de Vries equation, this solution is 
given by the Fredholm determinant 

w & = logdet(/ - tT 2 x) ); (1.16) 

see [10, p. 86, 174]. Tracy and Widom [28] introduced the cumulative distribution function 
F(x;t) = det(J - tT 2 {x) ) so that 

F(x;t) = exp(- [ (y - x)w(y;t) 2 dy); (1.17) 

J X 

in particular, F(x; 1) is the Tracy-Widom distribution. 

Operators for parts of the spectrum. In section 5 we show how Wi/ 3 and T( ) arise 
from the Airy group e %tT)Z on L 2 (R) where D = —i-^. By suitable changes of variable 
we arrange that the edge of the support of the equilibrium distribution is at zero and we 
consider the operators on L 2 (a, oo) that describe the probability that scaled eigenvalues 
lie in (a, oo). 

The relative positions of L 2 (R + ) and Wi/ 3 L 2 (H) are described by P + Wi/ 3 P + . Gen- 
erally we have Wt = e ttDi ' p_ e - ltD3 an d the complementary orthogonal projection is 

= I - W t = e itDZ 'P + e~ itDZ . For comparison, R + = T*P + T and R- = T*P-T 
are the Riesz projections on L 2 that have images H 2 and H 2 respectively. These formula? 
suggests an analogue of prediction theory such that the subspace T4 / t " L L 2 (R) corresponds 
to the subspace L 2 (R + ) which represents the future, and such that the unitary operator 
e itD 3 p} a y S the role corresponding to the inverse Fourier transform T* in the Hardy space 
theory of [17]. The projections r a W^T-a = e rfD3 P( Q , )0O )e~ ltD3 form a decreasing nest 
as a increases. The spectrum of the Hankel operator determines the limiting eigenvalue 
distribution via (1.13) and (1.15). 

The following table describes the analogy between the operators and subspaces in the 
various cases. 



G. Blower: Operators for soft and hard edges 6 



Classical Bulk 

Future Projection T*P + T T*P^- a ^ a) T 

Future space H 2 Ba/^L 2 

Subspace position e t2axjj2 c ^-2 

Painleve Equation cr-Py 

Hankel operator \l/ a 

Weyl Relations and Invariant Subspaces. 

Definition. Let (U t ) (t G R) be a Co (strongly continuous) group of unitary operators on 
an infinite-dimensional separable Hilbert space H, and let K be a closed linear subspace of 
H. We say that K is doubly invariant for (Ut) when U t K C K for all t G R. Further, K is 
simply invariant for (U t ) (t > 0) when U t K C if for t > and moreover n^ot/t-K" = {0}. 

Beurling and Lax characterized the subspaces of L 2 that are invariant for the shift 
operators S s : f(x) 1— > e lsx f(x); see [15, 17 page 114 ]. For notational simplicity, we 
sometimes write e %ix H 2 = {e %ix f(x) : f G H 2 }. 

A closed linear subspace T is simply invariant for the semigroup {S s : s > 0}, if and 
only if there exists a unimodular measurable function u such that T = uH 2 = {uf : / G 
H 2 }; such a « is uniquely determined up to a unimodular constant factor. In each case 
we start by making unitary transformations to identify u and to determine the relative 
positions of uH 2 and H 2 . Either uH 2 n H 2 = or there exist inner functions v and w, 
uniquely determined up to a unimodular constant factor, such that u = vw, uH 2 n H 2 = 
vH 2 and vwH 2 = vH 2 n wH 2 . In sections 4 and 5, we find wif 2 n -ff 2 = 0, so we factorize 
u(z) = E(z)/E(z) where E is a meromorphic function on C that has no zeros. By using de 
Branges's version of Beurling's theory [6], we are able to show that W from (1.1) is unitarily 
equivalent to r**r u » and hence that W is the reproducing kernel of some weighted Hardy 
spaces of holomorphic functions inside C+. 

Definition. A Weyl pair (U s , Vt) consists of a pair of Co unitary groups (U s ) (s G R) and 
(V t ) (t G R) on H that satisfy U s V t = e lst V t U s for all (s, t G R). 

The shifts S s (s G R) and the translations r t = e~ ltD (t G R) give a Weyl pair on 
L 2 ; moreover, this is the unique representation of the Weyl relations of multiplicity one on 
L 2 , up to unitary equivalence; see [32]. Katavolos and Power [16] obtained the following 
description of the invariant subspaces for a Weyl pair of multiplicity one. Let C be the 
space of orthogonal projections P onto closed linear subspaces K of 1? that are invariant 
under (S t ) (t > 0) and (r_ t ) (t > 0), where C has the strong operator topology. There is 
a homeomorphism p : {z : \z\ < 1} ^ £ such that : 



Hard Edge Soft Edge 

GtP+Ge e ^ 3 /3p +e -^ 3 /3 

Q1L 2 W^L 2 

ujl 2 n H 2 ^ e itx3 H 2 n H 2 = 

Pill Pn 

$1 r (0) 
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(1) p(-i) = 0, and p(i) = I; 

(2) p(z)L 2 with \z\ < 1 is simply invariant for both (S s ) (s > 0) and (r_ f ) (t > 0); 

(3) p{e l9 )L 2 with -tv/2 < 9 < tv/2 is simply invariant for (S s ) (s > 0) and doubly 
invariant for (r_ t ) (t e R); 

(4) p(e l0 )L 2 with — n/2 < n — 9 < n/2 is doubly invariant for (S s ) (s G R) and simply 
invariant for (r_t) (t > 0). 

In section 4 we introduce for the Jacobi ensemble an appropriate Weyl pair for the 
subspaces QiL 2 . On account of the natural ordering of the subspaces QtL 2 , and the 
probabilistic interpretation of (1.9), we naturally take the translations to be one of the 
groups in the Weyl pair; whereas we need to hunt down the other one. For the soft 
edge ensemble, the Weyl pair consists of the translations e lsD and the Schrodinger group 
e it{D 2 +x) ^ ag we dj scuss j n section 5. 

In section 6 we extend these ideas to a new context, namely the Mathieu functions. 
Here the KdV equation is 27r periodic and associated with an infinite-dimensional manifold. 
Whereas we do not propose that this corresponds to a natural random matrix ensemble, 
the results illustrate the scope of the theory of Tracy-Widom operators. 

2. Kernels from differential equations and the Marchenko integral equation 

In this section we prove some results concerning Tracy-Widom operators which are 
already known in specific cases from [8, 28, 29, 30]. Here B(if), c 2 and c 1 respectively 
denote the bounded, Hilbert-Schmidt and trace-class linear operators on Hilbert space H. 

Lemma 2.1. Suppose that A and B are bounded, measurable and real functions. Then 



W(x,y) 



A{x)B{y) - A(y)B(x) 
x-y 



(2.1) 



defines a self-adjoint and bounded linear operator on L 2 (R). 

Proof. This follows from the fact that Ma, M b and R + are bounded on L 2 . 



□ 



Proposition 2.2. Suppose that A and B are bounded, continuous and integrable functions 
such that 

d \A(x)~ 

dx~ B(x) = - 7 (a;) -a(x) B(x) ' K > 



where a, f3 and 7 are linear functions such that 

1 



C 



x 



y 



l{x)~l{y) a(x)-a(y) 
a(x)-a(y) (3(x) - (3(y) 





c 


a 




a 


b 



(2.3) 
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is a negative (semi-) definite constant matrix. 

(i) Then there exist continuous real functions F and G such that 



A(x)B(y)-A(y)B(x) 
x-y 



poo 

= / (F(x + t)F(t + y) + G(x + t)G(t + y))dt (2.4) 
Jo 



so P + WP + is a sum of two squares of Hankel operators. 

(ii) In particular when C has rank one, the operator P + WP + on L 2 (0, oo) is the square 
of a self-adjoint Hankel operator. 

Proof. We take a real bilinear pairing and write 



A(x)B(y) - A(y)B(x) = 



"0 -1" 




' A(x)' 




'MvY 


1 _ 




B(x) 




_B(y)_ 



(2.5) 



so by a short calculation 



\dx dy) 



A(y)B(x) _ i 


c a 




' A(x)' 




'MvY 


y \ 


a b 




_B(x)_ 







(2.6) 



We introduce a real symmetric matrix X such that X 2 = — C, let col [cos 9, sin 9] and 
col [— sin 9, cos 9] be the unit eigenvectors of X corresponding to eigenvalues Ai > and 
A2 > respectively. Then 



F(x) = X 1 (A(x)cos9 + B(x) sin 9), 
G(x) = X 2 (-A(x) sin9 + B(x) cos 9) 



(2.7) 



are bounded, continuous and integrable functions that satisfy 

(; 



d | d\ A(x)B(y)-A(y)B(x) 
.dy dxJ x — y 

= -F(x)F(y)-G(x)G(y) 



Vdx' + dxJ Jo { F ( x + t ) F (y + t ) + G ( x + t ) G (y + t )) dt - ( 2 - ? 



Hence both sides of (2.4) differ by f(x — y) for some differentiable function /; but both 
sides converge to zero as x or y tend to 00; so / = 0, and equality holds. 

□ 

This result enables us to calculate a determinant as in (1.4), (1.9) and (1.15). 
Theorem 2.3. Let A : (0, 00) — > R be a continuous function such that 
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J™uA(u) 2 du < 1. Then 

pOO 

W(u,v)= A(u + t)A(t + v)dt (2.9) 
Jo 

is the kernel of a trace-class operator on L 2 (0, oo) such that, when \k\ < 1, 

p OO 

K(x,z)-k 2 K(x,y)W(y,z)dy = kW(x,z) (2.10) 

J X 

has a solution K{x,z), which is a trace-class kernel, such that 
d 

— \ogdet(I-K 2 P {Xj0o) WP {Xj0o) ) = KK(x,x) (x>0). (2.11) 
In the subsequent proof, we shall use the self-adjoint Hankel operator r^) as in 

p oo 

T (x) f(t)= A(x + t + u)f(u)du (/Gl 2 (0,oo)), (2.12) 
Jo 

where T 2 x ^ on L 2 (0, oo) is unitarily equivalent to P( x ,oo)WP(x,<x) on L 2 (x, oo). 

Lemma 2.4. For \k\ 2 uA(u) 2 du < 1 there exists a solution L e L 2 ((0, oo) 2 ) to the 
integral equation 

rOO fOO 

L(x,s)~k 2 / L(x,y) / A(y + u)A(u + s) dudy = kA(x + s). (2.13) 

JO Jx 

Proof of Lemma 2.4. By the Hilbert-Schmidt theorem applied to rV \, there exist ((pj), 
an orthonormal basis of L 2 (0, oo), and real 7j such that A(x + y) = Yl^jLilj { Pj{ x ) { Pj{y) 
and 

°° poo poo poo 

V" 7 2 = / / A(u + v) 2 dudv= / uA(u) 2 du < oo. (2.14) 
• =1 Jo Jo Jo 

Our solution separates into series L(x,s) = Yl^LiXj( x )fj( s )^ where col[xj(x)] satisfies 
the equation with column vectors in £ 2 

[I - K 2 <S>(x)]co\[ Xj (x)} = KCo\[ im (x)} (2.15) 

and with matrix 

poo 

§(x)= -yjlk ¥j{y)¥k{y)dy . (2.16) 

L J M<j,k«x 
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By the Cauchy-Schwarz inequality, <&(x) defines a Hilbert-Schmidt operator on £ 2 with 
norm 

00 . roc poo 00 2 

j,fc=i Jx Jx i=i 

hence, 7 — k 2 $(x) is invertible whenever \k\ 2 \\$>(x) || c 2 < 1. We deduce that (2.15) has a 
unique solution and by orthogonality 

poo />oo />oo 00 

/ / \L(x,y)\ 2 dxds= / ^l^fcix 
Jo Jo Jo - =1 

/•OO 

00 _2 00 

J'=l J'=l 

□ 

Remark. When J2jLi < °°7 the operator is trace class and the determinant 

det(J — k 2 Q(x)) defines an entire function of k. 

Proof of Theorem 2.3. By Lemma 2.4, the integral kernels 

pOO 

W(x,y)= A(x + s)A(s + y)ds (2.19) 
Jo 



and 



poo 

K(x,z)= L(x,s)A(s + z)ds (2.20) 
Jo 

are trace class and satisfy the Marchenko integral equation 

P 00 

K(x,z)-k 2 K(x,s)W(s,z)ds = kW(x,z). (2.21) 

J x 

By iterated substitution, we deduce that 

K{x, z) = kW{x, z) + k 3 W(x, . )P (x , oo) {I - « 2 P(x,oo)WP( x , 00 ))- 1 P( x , 00 )W r ( . , z). 
We deduce, by simplifying the integral kernels, that 

kK{x,x) = k 2 ({I-k 2 Y 2 x) )- 1 A{x+ .),A(x+ .)> La(0>O o)- ( 2 - 22 ) 
By unitary equivalence as in (2.12), we have 

det(J - K 2 P iXi00) WP (Xi00) ) = det(J - k 2 Y 2 x) ) 
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and hence 

— logdet(/ - JP (x ,oo)WP {Xt0o) ) = - K 2 trace ((/ - ^r 2 ^)" 1 — r 2 x) ) , (2.23) 
where J^T^ is the rank-one operator that has kernel 

— V 2 (x) (u, v) = — I A(s + u + x)A(s + v + x) ds 
ox K ' ox Jq 

= -A(u + x)A(v + x); (2.24) 

hence 

^ logdet(J - k 2 T 2 x) ) = -k 2 {{I - K^l^Aix + .), A(x + .)> La(0>oo) - (2-25) 
By comparing the terms in the power series in k, we deduce that 

kK(x,x) = ^logdet(I - k 2 P {Xj00) WP {Xj00) ). (2.26) 

□ 

Corollary 2.5. Suppose further that A is an entire function such that J °° u\A(z+u)\ 2 du < 
oo for each z G C, and let T( 2 ) be the Hankel operator on L 2 (0, oo) that has kernel 
A(z + s + t). Then 

^logdet(/- K 2 r 2 2) ) (2.27) 
defines a meromorphic function on C, extending kK(x,x). 

Proof. By Morera's theorem, z h- > defines an entire function with values in c 2 , and 
hence det(J — k 2 T 2 z ^ defines an entire function. The formula (2.27) defines a holomorphic 
function, except at those isolated points where the determinant vanishes, and these give 
rise to poles. 

□ 

In some cases y = -^K{x,x) satisfies a Painleve equation as in [14, p 344]; that is, 
y" = F(y', y, x) where F is rational in y and y', and analytic in x, and such that the only 
movable singularities of y in C are poles. In particular, the bulk kernel gives rise to the 
o form of Py, the hard edge ensemble gives rise to the Pjn equation and the soft edge to 
P n as in [28, 29, 30, 12]. 

The Painleve ODE test asserts that every ordinary differential equation that arises 
from a partial differential equation via a Marchenko linear integral equation may be trans- 
formed to a Painleve equation; see [1]. 
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A special feature of (2.10) is that W is the square of a self-adjoint Hankel operator V 
on L 2 (0, oo), and Proposition 2.2 gives an explicit construction of the symbol for T. Up to 
unitary equivalence, this holds under some general spectral conditions which we list below. 
In section 3 we give a means for verifying (hi) for operators of the form (1.1). 

Proposition 2.6. Suppose that W is a linear operator on H such that 

(i) the nullspace of W is either trivial or infinite-dimensional; 

(ii) W is not invertible; 

(in) W is bounded and self-adjoint, and W > 0; 

(iv) W has a simple discrete spectrum. 
Then there exists a self-adjoint Hankel operator T on L 2 (0, oo) and a unitary operator 
U :H -> L 2 (0, oo) such that W = U*T 2 U. 

Proof. Megretskii, Peller and Treil [21, p 257] have obtained sufficient conditions for a 
self-adjoint operator to be unitarily equivalent to the modulus of a Hankel operator. Under 
the more stringent condition (iv), their construction gives a Hilbert space K, a bounded 
linear operator X : K — > K, and vectors £,rj G K such that the Hankel operator 

Tf(t)= h(s + t)f(s)ds (/GL 2 (0,oo)) (2.29) 
Jo 

with symbol h(t) = {e tx ^,rf) K is unitarily equivalent to IU 1 / 2 ; thus V is realized from a 
balanced linear system in continuous time with one-dimensional input and output spaces. 

3. Reproducing kernels and the bulk of the spectrum 

In this section we recover the bulk kernel as a reproducing kernel, and show more 
generally why operators of the form (1.1) are positive on weighted Hardy spaces. Let E 
be a meromorphic and zero-free function on C and let E*(z) = E(z), which has similar 
properties. We also introduce the meromorphic functions A(z) = (E(z) + E*(z))/2 and 
B{z) = (E*(z) — E(z))/(2i), which have A(x) and B(x) real for real x. 

Let EH 2 be the weighted Hardy space of meromorphic functions g on C_|_ such that 
g/E belongs to the usual Hardy space H 2 , and with the inner product 

r°° _ dt 

{9u92)Em = (gi/E,g 2 /E) H 2 = J gi(t)g 2 (t) jj^yp - t 3 - 1 ) 

Similarly we can introduce E*H 2 . When £ e C+ is not a pole of E, the linear functional 
g i— > g(Q is bounded on EH 2 , and hence given by g(() = (g, k^) eh 2 ? where the reproducing 
kernel is 

= mm 

° 2m(C-z) ' 
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We introduce O as the domain consisting of points z G C+, that are not poles of E or 
E*. Let u(z) = E*(z)/E(z), which is meromorphic and unimodular on the real line, let 
M u : EH 2 -> E*H 2 be the isometry M u f = uf, and let T s : H 2 -> H 2 be the Toeplitz 
operator T u = R + M U R+. 

Theorem 3.1. (i) The operator W on EH 2 that has kernel 



, E*(z)E*(w)-E(z)E(w) . tnn . 

w(z ' w) = 2„ ( ;- ffl ) <*■ " e n) (3 - 3) 

compresses to an operator EH 2 — > i?i7 2 that is unitarily equivalent to r^T^, where Fu : 
H 2 — ► H 2 is the Hankel operator F^ = R- M u * R+ . 

(ii) There exists a unique Hilbert space H(W) of holomorphic functions on O such 
that W(z, w) is the reproducing kernel for H(W). 

(in) Suppose that T a has a non-zero nullspace K. Then F n restricts to an isometry 
K — > If 2 . 

Proof, (i) We write 



/: 



E*(z)E(t)-E(z)E(t) f(t)dt E(z) f°° f(t)/E(t) ^ 



2m(z-t) E(t)E(t) 2ni J_ 00 t-z 

E*(z) r u*(t)f(t)/E(t) 



2ni 

and hence by Cauchy's integral formula we have 



f (3.4) 



Wf(z) = E(f/E - M U R+M U * (f/E)) 

= EM u R.M u *(f/E). (3.5) 

The map V : EH 2 — > H 2 : f i— > f/E is a unitary equivalence with adjoint V* : g \- > 
Eg, and T^r s : H 2 ^ H 2 reduces to = R + M U R_M U *R + , so (Wf,g) EH 2 = 

(V*FlF u Vf,g) H2 for all /, <7 £ EH 2 . 

(ii) By (i), is a positive operator on EH 2 , so the kernel W(z, w) is of positive type 
on Q; further, z i— > w) and w i— > w) are holomorphic on O. Hence we can apply 
[2, Theorem 2.3.5] to obtain the Hilbert space of holomorphic functions such that W(z, w) 
is the reproducing kernel. 

(iii) By [23, p 89], we have T?T a = f - F* u F a , which leads directly to the identity 
K = {f EH 2 -.1^11 = \\f\\}. 

□ 
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Corollary 3.2. Suppose that u belongs to H°° so that E*H 2 is a closed linear subspace 
of EH 2 , and let K = EH 2 E*H 2 be the orthogonal complement of the range of M u : 
EH 2 — > EH 2 . Then K equals H(W) and has reproducing kernel 

k{x) = M*)b{»)-bmm») m 

7T{W — Zj 

Proof. First, one can check by calculation that 



, s ttt/ v E*(z)E*(w)-E(z)E(w) 
K w (z) = W(z, w) = I (3.7) 

2ni(z — w) 



Then we observe that 



2m(z-w) v ' 2m(z-w) 
lies in the range of M u ; so for g E K the proof of Theorem 3.1(i) simplifies to give 

(g,K w ) EH2 = (g,k w ) EH2 = g(w) (weO). 

□ 

Bulk of the spectrum. Thus when u is an inner function we can identify H[W) 
explicitly as the orthogonal complement of a shift- invariant subspace of EH 2 . In particular, 
by taking the entire function E{z) = e~ taz , we find u{z) = e 2iaz and the reproducing kernel 
for K = EH 2 e E*H 2 to be 

TT(Z - W) 

as in the bulk kernel B a / 7T (z 1 w) of (1.3). Here we have EH 2 = T*I?\— a, oo), and \l/ a = 
T*\L 2 \— a, a] gives a unitary isomorphism L 2 [— a, a] — > K with ^/a^a = ^o/tt- The Hankel 
operator T s is isometric on H 2 Q e 2iax H 2 ~ if. 

The Paley-Wiener theorem [17, p. 179] characterizes K as the space of functions 
/ G L 2 (R) that are entire and of exponential type with 

lim sup log \f(iy)\ < a. (3.10) 

Alternatively, we can characterize the subspaces by their scaling properties. Let (5 t ) (t £ 
R) be the unitary dilatation group on L 2 (R) with S t f(x) = e*/ 2 /(e^). In [16], Katavolos 
and Power characterize the lattice of closed linear subspaces of L 2 that are simply invariant 
for both S s (s > 0) and 5 S (s > 0). 
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Proposition 3.3. The closed linear subspace B t L 2 



is simply invariant for (S s ) (s < 0), 



doubly invariant for (r s ) (s G R) and invariant under J. Conversely, if K is any closed 
linear subspace of L 2 that is simply invariant for (5 t ) (t < 0), doubly invariant for (r s ) 
(s G R) and invariant under J, then K = B a L 2 for some a > 0. 

Proof. We have S- s = JF*5 S JF and r s = JF*5'_ S JF, so we shall characterize the subspaces 
L 2 [— 7rt, lit] under the operation of S s , S s and J. Now L 2 [— tw, irt] is clearly doubly invariant 
for (S s ) (s G R), and 8 S L 2 [— nt, irt] = L 2 [—nte~ s , ixte~ s ]; so L 2 [— tn, irt] is simply invariant 
for (5 S ) (s > 0). Conversely, all closed linear subspaces K of L 2 that are simply invariant 
under (5 S ) (s > 0) and doubly invariant under (S s ) (s G R) have the form K = L 2 (— a,6) 
for some a, b G R U {oo} by a simple case of Beurling's theorem. When K is additionally 
invariant under J, we need to have a = b; hence K = L 2 [— a, a]. 

4. Hard-edge operators and Sonine spaces 

In this section we consider the operators for the hard edge case and associated sub- 
spaces. Let J u be the Bessel function of the first kind for real v > —1/2, and let 



which is entire and of order 1/2 as in [11, p 190]. Then E(z) = l/h(z) is a meromorphic 
function, with no zeros, such that 



We recognise the first factor on the right-hand side from (1.8), and the left-hand side from 
(4.4); but Corollary 3.2 does not apply directly to E*(z)/E(z); so we introduce operators 
that correspond to these kernels indirectly by means of the Hankel transform as in [25, p 
298]. The Hankel transform of / G L 2 (xdx; (0, oo)) is 



and the unitary operator U : L 2 {xdx, (0, oo)) — > L 2 (R) by Ug(^) = e ^g(e ^) such that 



E*(z)E*(w) - E(z)E(w) 
2iri(z — w) 




). (4.2) 




(4.3) 



On L 2 (xdx; (0, oo)) we introduce the unitary dilatation group (5 t ) by 5 t g(x) — e t g(e t x) 



U*r t U = 5 t . 



G. Blower: Operators for soft and hard edges 16 



Lemma 4.1. Let Gi be the integral operator on L 2 (R) that has kernel function 

e-^-V^e-'-*-"). (4.4) 

Then Gt is a self-adjoint and unitary operator such that G\ = I, and G(jt = T- t Gg. 

Proof. From the shape of the integral kernel, the identity G(U = T-gUl-L v is evident. 
Further, Hankel's inversion formula leads to the identity 7i 2 = I, whence to 

G e UU*G e = T- t UH v H v U*T t = I. (4.5) 

The identity (4.4) is evident from the definitions, and by (4.5) is equivalent to the 
scaling property H. v 8t = 8-{H. v of the Hankel transform as in [25, p. 299]. 

□ 

The following result on position of subspaces contrasts with Corollary 3.2. Here V 
denotes Euler's gamma function. 

Theorem 4.2. (i) The operator Qi = G^P+Gi on L 2 (R) is an orthogonal projection, 
(ii) The range of TQtT* equals e lS - x u v H 2 , where the meromorphic function 

is holomorphic on {z : < 0}, and unimodular and continuous on R. 

(in) Whereas u* v H 2 n H 2 = {0}, for v > the subspace K = (uJT 2 ) n H 2 is non-zero, 
and : H 2 — » H 2 restricts to an isometry K — > H 2 . 

Proof, (i) This follows directly from the Lemma. 

(ii) Our aim is to show that the range of the orthogonal projection T*QiT is simply 
invariant under the S\ for A > 0. By Plancherel's theorem we have 

S X FQ £ L 2 = fr- X G e P + L 2 = FG T X+e P+L 2 , (4.7) 

where T\ +£ P + L 2 = L 2 (A + oo) C L 2 (£, oo) and n\ >0 L 2 (\, oo) = 0. Consequently 
by Beurling's theorem, there exists a unimodular and measurable function u v such that 
TQqu v L 2 = u u H 2 , and u v is unique up to a unimodular constant factor. One can easily 
deduce that TQtL 2 = e lix u v H 2 . 

The Fourier conjugate of Qi is TQtT* = TGiT*TP+T*TGiT*, wherein we recognise 
FP + F* as R- : L 2 — > H 2 . To determine the range of TQ^T* , or equivalently the subspace 
J r G^L 2 (0, oo), we write 

TG t f{x) = ^ e- ix * [°° e^-*-" Me-'-t-^fWdrj ^ 

J -oo J0 
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for / G L 2 (0, oo), and then reduce this integral by simple transformations to 

/oo 
e -(l+u+ix^ e u^j^ e -^ d ^ (48) 
-oo 

The substitution y = reduces the final integral in (4.8) to a standard Mellin transform 
[25, p. 263], and we identify u v from 

(hi) Let E v (z) = e~ iz log ^r((l + v - iz)/2) so that E v is meromorphic and zero-free 
with simple poles at — i — vi — 2ki for k = 0, 1, . . ., and u v {z) = E*(z)/ E u (z) has simple 
zeros at zt = —i — vi — 2ki for k = 0, 1, . . . and simple poles at i + vi + 2ki for k = 0, 1, . . .. 
The function u u (z) is holomorphic in the lower half plane, but does not define a bounded 
holomorphic function on {z : Qz < 0} since the series Ylh=o^ z k/(^ + \ z k\ 2 ) diverges, 
violating Blaschke's condition for the zeros of a non-trivial function in H°° or H 2 as in 
[17, p. 92]. Hence the equations hi{z) = ul(z)h 2 (z) with hi, h 2 G H 2 has only the trivial 
solution hi = h 2 = 0; so u* v H 2 n H 2 = 0. Note that log |it£(z)| is subharmonic on the C+, 
but is not the Poisson integral of a measure on R. 

We take a > and v + 1/2 > A > 1/2, and let 

/(x) = a-^x 1 ^-!^ - a 2 )^" 1 )/ 2 ^! («v^)l (aiCo) (x), 
with Hankel transform 

g(t)=t 1 / 2 H l/ (x- 1 / 2 f(x);t). 

Then by a result of Sonine [6 p. 301 , 24 p. 75, 26 p. 38], both / and g are supported on 
(a, oo) and we have 

/>oo /»oo 

/ g(t)t-^ 2+ix dt = u v {x) I /(£)£- 1/2 - ia: eft (ieR). (4.10) 

Hence when a = 1 there exist non-zero functions hi,h 2 G -ff 2 such that h 2 (x) = u u {x)h\{x), 
so /12 £ u u H 2 . Now we apply Theorem 3.1(iii) to deduce that Y Uv \H 2 nu u H 2 is an isometry. 

□ 

Proposition 4.3. (i) The Hankel operator 5v = P+GtP + on L 2 (0, oo) has $ 2 = P + QgP + . 
(ii) The operator &g on L 2 (0, oo) is Hilbert-Schmidt, and each non-zero 
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/ G L 2 (xdx, (0, 1)) such that 



Jo 



dy 



(4.11) 



(4.12) 



corresponds to an eigenfunction g G L 2 (0, oo) of <Ev with eigenvaiue ^Ay^- 
(iz'i) The kernel of Qi as an integral operator on L 2 (R) is 

e-'-ZMe-t-^e-u-^JKe- 1 -**) - e-^-^Jlje-^e-^Me-^) 

e -2t-2i _ e -2£-2r, 

(iv) det(I - zF°> a ) = det(I - z$ 2 a) ) for a = -(1/2) loga and a > 0. 

Proof, (i) For t > we have the Hankel condition <£>£T t = r t *$^, where here (r t )t>o denotes 
the semigroup of translation operators on L 2 (0, oo). Then one uses Theorem 4.2(i). 

(ii) The kernel function is clearly symmetric, real-valued and square integrable, since 



pco poo pC 

/ / e -2(^) Jv ( e -(^+0)2^ = 
Jo Jo Jo 



ue 



-2i-2u 



Ju(e~ e ~ u ) 2 du < oo (4.13) 



due to the asymptotic formula J u (x) x x u /T(u+l) as x — > 0+. Hence gives a self-adjoint 
operator of Hilbert-Schmidt type. The operator U restricts to a unitary L 2 {xdx; (0, 1)) — > 
L 2 (0, oo), and under this transformation the eigenfunction equations correspond via #(£) = 
e-«/(e- 2 «). 

(iii) We use the method of proof of Proposition 2.2 to verify the stated formula for 
Qi = GgP + Gi, which is essentially the square of a self-adjoint Hankel operator. With 
A(0 = e-*J v (e-Z) and = e" 2 ^ J^e^), we have 



d_ 



where 



-2? 



hence 



d_ 



dt]J 



-2? _ 



A' 






B 




> 


-1 




-1 


-v 2 ) 


-1 


e 


2»? _ 


e~ 











- e 


2^7 





-1 





"A" 







+ 





-1 



-1 (e" 2 " - v 2 ) 
-1 -1 



(4.14) 



-1 

1 



+ 



(4.15) 



8_ d_ 
di drjJ J 



A(£ + u)A(rj + u) du. 



Hence we can obtain the stated identity by following the proof of Proposition 2.2. Alter- 
natively, one can transform a formula in [25, p. 303]. 
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(iv) The unitary equivalence between T 2 ((0, 1), dx) and L 2 (0, oo) involves g{x) i— > 
\/2e~£g(e~ 2 ^), so i^ ' 1 ) is unitarily equivalent to the operator that has kernel 

2e _^ F(0 , 1)(e _ 2 ^ e _ 27;) = e-ZJAe-t)e- 2 nJl(e-n) - e^je^e^ J v ( e -") ^ 

e -2£ _ e -2ri 

which we recognise as the kernel of & 2 y Comparing the spectra of the compressions to 
L 2 (0, a) and T 2 (a, oo), we deduce that 

det(J - zF°' a ) = det(I - zP {a ^\ 0) P {ct)Oo) ) = det(J - z$ (0 )P( a ,oo)$(o))- (4-16) 

Finally, $(o)-P(a,oo)^ > (o) equals ® 2 a ^ since they both have kernel 

poo 

/ e- ? - u J ! ,(e- ? - u )e- ?7 - u J ! ,(e- ?7 -")rfw. (4.17) 

J a 

□ 

Theorem 4.4. Let T be the operator 

= -| (e 2 ^) + (- 2 " l)e 2 V(0- (4-18) 

(i) Then T is an essentially self-adjoint and positive operator on C^°(R) in L 2 (R), so 
that V t = e - M T- u / 2 (t G R) deques a C group of unitary operators on L 2 (R). 

(ii) The unitary groups (V^) s€ r and (r t ) te R satisfy y s r t = e* s V t F s for s,t G R. 

(iii) The subspace QiL 2 is doubly invariant for (V s ) with s6R and simply invariant 
for (r_ t ) for £ > 0. Conversely, if if is a non-trivial closed linear subspace of L 2 that is 
simply invariant for T-t (t > 0) and doubly invariant for V s (s G R), then K = Q a L 2 for 
some real a. 

Proof, (i) The simplest way of proving that the operator T is self-adjoint is to compute 
its spectral resolution. By simple transformations of the Bessel equation [14, p 171], we 
have 

-e 2 «(|^+2^ + ^-l)(e^-^^ (4.19) 

so that e~ ? ~^~ T ' J I/ (e~^~ i ~' n ) is an eigenfunction of T corresponding to the eigenvalue 
e -2e-2r] > o HankePs inversion theorem [25, p 299], the functions XyJ iy (\xy) give 
a complete spectral family in L 2 (xdx; (0, oo)), and the unitary transformation U takes 
XyJ v {Xxy) to e~ ? ~^~ T ' J I/ (e~^~ i ~ v ) after an obvious change of variable. By Stone's theo- 
rem, (— i/2) logT generates a C unitary group T~ ls l 2 . 
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(ii) We have 

/CO 
e-i-v-tJvie-t-t-^fWdn 
-co 

/CO 
e i871 e-^- £ - v J v {e-^- e - v )f{r]) drj 
-co 

= G e SJ(0; (4.20) 

hence G(V s Gi = S s . When we conjugate the familiar Weyl-von Neumann relation T-tS s = 
e lst S s r- t by G e ® G\ we obtain G e T- t G e G e S s G e = e ist G e S s G e G e T- t G e or r t V s = e lst V s r t . 

(iii) From earlier relations, we have 

V S Q £ = V s G e P+G e = G £ S S P+ = G e P+S s = G e P+G e G e S s = Q e G e S s , (4.21) 

which shows that the range of Qe is mapped onto itself by V s ; further 

T-tQe = T-tGgP+Gg = GiTtP+Gg = G £ P[ tt00 )T t Ge = G^P[ t)00 )G^r_ t , (4.22) 

has range contained in the range of Qt for t > 0, so QtL 2 is simply invariant. 

To obtain the converse, we consider the Fourier transforms of the groups. On 
e lix u v Pi 2 , the unitary semigroups operate as 

V s =TV S T* :f(x) ^e iis u l/ (x)u l/ (s-x)f(x-s) (sgR); (4.23) 

Tr-tT* = St- To verify (4.23), we recall the flip map J by Jf(x) = f(—x), and observe 
that TT = J, and J 7 * J 7 * = J. We have 

TV S T* = TGzSsGzT* = SiM u V T*S S TT* GiT* (4.24) 

so that TV S T* = S^M Uv t s JSiM Uv J. Using the Weyl-von Neumann relation for r s and 
Sg, one can easily simplify this expression to obtain V s = TV S T* = e lis M Uv N s t s , where 
N s f(x) = u v (s — x)f(x). The functions u v satisfy u v {— x)u u (x) = 1 and 

e u„[x)u u {s x)-e Z r{{1 + u _ lx) /2)r{{1 + u + lx _ ls) /2 y 

Suppose that K is such an invariant subspace. Then by Beurling's theorem, there 
exists a unimodular and measurable function w such that TK = wH 2 ; further, this w 
is uniquely determined up to a unimodular constant multiple. We apply TV a T* to this 
identity, and deduce by double invariance and (4.23) that 

{e Us u v {x)u„{s - x)w(x - s)f(x - s) : f e H 2 } = wH 2 ; (4.26) 
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so that, 

e iis u u (x)u u (s — x)w(x — s) = c(s)w(x) (s G R) 

holds for some c(s). We re-arrange this to u v {s — x)w(x — s) = u v (—x)w(x)e s c(s), 
then solve to obtain w{x) = e lolx u v {x) for some a G R. Hence TK = e iax u u (x)H 2 , so 
K = Q a L 2 by Theorem 4.2(h). 

□ 

Corollary 4.5. The unitary groups (V s ) (s G R) and (S t ) (t G R) form a Weyl pair. The 
closed linear subspaces that are simply invariant for V s (s > 0) and S t (t > 0) have the 
form 

K a ,p = {e iax ~ il3x \ v {x)f{x) :fEH 2 } (a G R,/3 > 0). (4.27) 



Proof. We shall check simple invariance of K a ^ under V s (s G R), since the other 
statements follow easily from the proof of Theorem 4.4 and the Katavolos-Power theorem. 
We have 

V a : e iax - llix \ v {x)f{x) ^ e Us - ias -^ s2 e tax -^ x2 u u (x)e 2 ^ sx f(x - s), (4.28) 

where e 2lf3sx f(x - s) G S 2f 3 s H 2 C H 2 for 8 ,0 > 0. 
5 Soft edge operators and the Airy group 

With D = —i-§^ the Airy group e ltD3 is a Co group of unitary operators, as defined by 

e itD3 f(x) = -= / e^+* x Tf{Z) df. (5.1) 

In this section we shall consider how the Airy group is related to the kernel function of 
(1.11). Here J t denotes the operator e ltL)3 J on L 2 (R), not a Bessel function, and we shall 
use a subscript t to indicate scaling of the space variables x and y with respect to time t. 

Lemma 5.1. The operator J t = e ltI)Z J is self-adjoint with J 2 = I, and J t as an integral 
operator on £ 2 (R) has kernel 

1 / x + y \ 

w 7 ^" (5 ' 2) 

Proof. For a compactly supported and smooth function / we have 

Je* tD3 Jf(x) = -= / e'^-^/CO d£ = e-* tD3 f(x), (5.3) 
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so J 2 = I. Further, the kernel of J t is given by 

-1 />oo />oo 

"Jf(x) = - < ■ / e^mdydt 

^ J — oo J — oo 



AtD 6 



-oo 

"OO 



JJ_ / Jft+ii(x+y) 



oo k Z7r </-oo ' 



Since the Airy function is real valued, it also follows that J t is self-adjoint. 

□ 

Most of the next result is essentially contained in [28, Lemma 2], but we include a 
proof for completeness. 

Proposition 5.2. (i) The operator 

W t = e UD3 P-e- ttDi = J t P + J t (5.5) 

on L 2 (R) is an orthogonal projection and the range ofTWtT* equals e w * H 2 . 
(ii) The Hankel operator Fqj = P + JtP+ has square Tq t = P + WtP+. 
(in) The kernel ofW t as an integral operator on L 2 (R) is 

W ^ V) = J^TTs Wi/s ( (3^73 , j^tjs) (5-6) 
where Wi/ 3 is the Airy kernel as in (1.11). 

Proof, (i) By Lemma 5.1, we have W 2 = J t P + J 2 P+Jt = JtP+Jt = W t , so that W t is a 
projection; further W£ = W t . The range of W t equals {W t k : k G L 2 (R)}, or equivalently 
the range of JtP+- 

If / G L 2 (R+), then Ff(£) = (5(0, where G G H 2 . Since e~ itD3 is unitary, we have 
W t L 2 = e ltI)Z P_e~ ltc>3 L 2 = e ttL)3 P-L 2 , and hence the image of WL 2 under the Fourier 
transform T is TW t L 2 = {e*** 3 F(£) : F G H 2 }. 

(ii) We have r ,t = P + e lt ° 3 JP + and hence 



T 2 it = P + e UD JP+e UD JP+ 



= P + e ttD3 JP+JJe ltD3 JP + 

= P+e itD3 P-e~ ltD3 P + = P+W t P+. (5.7) 

(iii) It also follows from the Lemma that the kernel function is 

w '^ )= w*L Ai ((3^) Ai ^)*' (5 - 8) 
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a formula which reduces to (5.6) and (5.7) on account of the identity 

Ai(x)Ai'(y) - Ai'(a;)Ai(y) 



poo 

/ Ai( 

JO 



x + u)Ai(u + y) du = 



x-y 



(5.9) 



This formula is presented by Tracy and Widom in [28], and follows from Proposition 2.2 
since A(x) = Ai(x) satisfies 



d 


'A' 




"o r 




A~ 


dx 


B 




x 




B 



□ 

Evidently Wq = P_, and the relative positions of the ranges of P- and Wt are de- 
scribed in the following Proposition, with different conclusions from Theorem 4.2. 
Definition. [9] A function G G H 2 is said to be cyclic (for the backward shifts) when 
span-fiS^C; t > 0} is dense in H 2 . Likewise, / G L 2 (R+) is cyclic when span{r t */ : t > 0} 
is dense in L 2 (R + ); g G L 2 (R_) is cyclic when span{r t *(7 : t < 0} is dense in L 2 (R_). 

Proposition 5.3. (i) For each t ^ 0, the subspaces W t L 2 nL 2 (~R_) and (W t L 2 ) ± nL 2 (R + ) 
equal {0}; while any non-zero vector in W t L 2 n L 2 (R + ) or (W t L 2 ) ± fl L 2 (R_) is cyclic, 
(ii) For each t > 0, the operator W t on L 2 (R_) © L 2 (R + ) has block matrix form 



P-W t P- P-W t P+ 
P+WtP- P+W t P+ 





"B c 2 ' 


G 






c 2 c\ 



(5.10) 



(in) For any real t, the operators P + W t P- and P-W t P + are Hilbert-Schmidt. 

Proof, (i) First we check that W t L 2 fl L 2 (R_) = 0, or equivalently by Proposition 5.2(i) 
that e u ^H 2 n H 2 = 0. Suppose that F,G G H 2 are non-zero and satisfy e^ 3 P(0 = G(f ) 
for almost all ^ G R. Then K(() = e lt < F(Q - G(() is a holomorphic function with zero 
boundary values at almost all points of R; so by the Lusin-Privalov theorem, K(() is 
identically zero on C + . Now by Szego's Theorem [17, page 108], the integrals 



\og\F(£ + in)\ 



d£ and 



log | G(f + i77)| 



converge. But this contradicts the identity e 1 ^ F(() = G((), since 

r0 ° ^ + tv) 3 



(5.11) 



i+e s 



-de 



(5.12) 



diverges for 77, t > 0; so P = G = 0. Likewise the only solution of the equation e lt ^ P(£) 
G(0 with F,GeH 2 is F = G = 0. 
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Next prove that all non-zero vectors in W t L 2 flL 2 (R + ) are cyclic. Suppose that G^O 
is a non-cyclic vector in H 2 fie**? 3 H 2 so that G(£) = e lt ^ 3 F(£) for some F e i? 2 , and where 
G _L wi7 2 for some inner function u. We have G H 2 ; so we introduce inner functions v 
and to, and an outer function 6>, such that uG = v6 and F = w9. Then, as in [9, Theorem 
3.1.1], 

e ue = % = — (5-14) 
F uw v ' 

is a quotient of inner functions and hence is of finite Nevanlinna type, but the corresponding 

logarithmic integral (5.12) diverges, and we have a contradiction. (The author conjectures 

that W t L 2 fl L 2 (R + ) = so that W t L 2 and L 2 (R + ) are in general position, since any 

non-zero elements in the intersecting subspaces would satisfy some implausible equations.) 

(ii) The Hankel operator r 0) t = P+J t P+ = P+e ltD * JP + has kernel 

-_I (0>TO) (x)Ai(^t^)l (0iOO )(y), (5.15) 

which is of Hilbert-Schmidt type; see [23, page 46]. Indeed, the integral 



o 



1 / x + y 



,-Ai 



dxdy 



(3t) 2 / 3 \{3ty 

may be transformed by the substitution u = x + y to the convergent integral 

1 f°° ( u \ 2 

m^L wAi w^ du<oo; (5 ' 16) 

here we use the bounds from [11, page 43] 

Ai(x) = ^171(1 + 0(x-" 2 )) exp(-^ 3 / 2 ) (x - 00). (5.17) 

Hence the off-diagonal operators P-WtP+ = P-Jt(P+JtP+) and 
P + WtP- = (P+JtP+)JtP- are Hilbert-Schmidt. For the bottom-right entry, we have a 
stronger conclusion, namely that P + WtP+ = (P+JtP+)(P+JtP+) is trace class. 

(iii) When we replace t > by t < 0, we need to switch the roles of P+ and P- in 
the previous discussion and we deduce that P-W t P+ and P + W t P- are Hilbert-Schmidt, 
while P-WtP- is of trace class. 

□ 

Theorem 5.4. (i) The C unitary groups S s and U t = e ~ lt ( D ~ x2 ) satisfy the Weyl relations 
S s U t = e lst U t S s for s,teR. 
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(ii) For a > and real 5, the subspace e lx3 / 3 - tax2 + t5x H 2 j s simply invariant for S s 
(s > 0) and U t (t > 0). Conversely, if T is a non-trivial simply invariant subspace for S s 
(s > 0) and for U s (s > 0), then T = e ix 3 /s-iax 2 +iSx H 2 for SQme a > and real 5. 

Proof, (i) One can prove directly that the operators U t defined by 

U t f(x) = e t(x2t - xt2+t3 ^f(x - t) (5.18) 

define a Cq unitary group on L 2 (R). This formula for the Ut of (5.18) was obtained by 
the method of characteristics. Indeed, when / is differentiable, the function g(x, t) = 
e i(x 3 -(x-t) 3 )/3j^ x _^ sa ti s fles 

dg dg . 2 
ot ox 

g(x,0) = f(x); (5.19) 

and so Utf(x) = e - lt ( D ~ x ^f(x) = g(x,t) gives the unique solution of the initial value 
problem (5.19). 

Let V be the unitary operator V : f(x) i— > e txS / 3 f(x) on L 2 (R), then clearly S s = 
V*S S V. The generator of the unitary group V*UtV equals 

-iV*(D - x 2 )V = -i e - tx3 / 3 (D - x 2 )e lx3 / 3 = ~ = -iD; (5.20) 

ox 

so by the uniqueness of groups with given generator we have V*U S V = e~ lsD = t s and 
hence U s = Ve~ lsD V* = Vt s V*. By conjugating the Weyl relations r s S t = e~ ist S t r s for 
(s,t £ R) by V, we can deduce (5.19). 

(ii) Clearly any T = e lx3 /^- %ax2 + tSx H 2 is simply invariant under S a (s > 0), and we 
can use the preceding calculations to show that T is also invariant for U s (s > 0). Indeed, 
for g G T we can take / G H 2 such that g(x) = e * x3 /3-tax 2 +i8x anc j we 

U s g = u s (e lx3/s - iax2+l5x f) 
= Vr s V*V{e- iax2+iSx f} 
= Vr s {e- iax2+lSx f} 

g2iasx— icts 2 — iSs ^ix 3 /3— icex 2 +i5x j ^ g^j 21) 

where f(x — s) is an H 2 function, so U s g belongs to the subspace S2 as T of T. This proves 
the forward implication. 

To prove the converse, we take any T that is simply invariant as in the Theorem, and 
observe that V*T is simply invariant under S s (s > 0) since V* commutes with S s , and 
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V*T is also invariant under r s (s > 0) since t s V*T = V*U S T C V*T. By the Katavolos- 
Power Theorem [15], there exist a > and a real 6 such that V*T = e ~ tax2 +lSx H 2 , and 
hence T has the required form. 

□ 

Corollary 5.5. The unitary groups r_ s = e tsD and U t = e **( D2+x ) form a Weyl pair on 
L 2 , such that the lattice of subspaces {T-sU- a Wi/ 3 L 2 : S G R, a > 0} gives the set of 
simply invariant closed linear subspaces for (r_ s ) (s > 0) and (U t ) (t > 0). 

Proof. The strongly continuous unitary semigroups r_ s = ^SgJ 7 (s > 0) and J-*UtT 
and T* StT (t > 0) have generators iD and z(D 2 + x) respectively, so (r_ s ,?7t) forms a 
Weyl pair by Theorem 5.4(h) and Proposition 5.2(i). Under the Fourier transform, we 
have 

^r_ < 5r_ a2 f7_ Q W 1/3 L 2 = SsS^U^VH 2 

hence the subspaces correspond as stated in the Corollary. 

6. Mathieu functions and periodic potentials 

One can construct Tracy-Widom operators over the circle by means of the differential 
equations of section 2. Kernels of the form (6.6) below arise in the theory of random 
unitary matrices with Haar probability measure as in [22, p. 195]. The purpose of this 
section is to introduce examples beyond the list in [30]. 

Let S be the 2x2 fundamental solution matrix of Hill's equation with smooth n- 
periodic potential g, so that 



then det S = 1, and A (A) = trace S(n) defines the discriminant. When A is real, evidently 
A(A) 2 > 4 if and only if the eigenvalues of S(n) are real, and A(A) 2 = 4 occurs if and only 
if S is periodic with period it or 2%. The periodic spectrum 

A = {A < Ai < A 2 < A 3 < A 4 < . . . < A n / oo} (6.2) 

of Hill's equation consists of those real A such that 



dx 



1 

-(X + q(x)) 



S(0) = 



y" + (X + q)y = 



(6.3) 
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has a non-trivial ix or 2tt periodic solution as in [18, p. 11]. The discriminant satisfies 

4 - A 2 (A) = 4(A - A ) A ^- 1 ~ X l^- X \ (6.4) 

3 = 1 3 

Theorem 6.1 For each real a there exists an infinite sequence of X n such that Hill's 
equation (6.3) with potential a cos 2a: has a non-trivial 2-K-periodic and real solution A. 
For such A, let W be the kernel 

sm(x — y) 

which continuously differentiable and doubly periodic with period 2tv. 

(i) Then W defines a self-adjoint and Hilbert-Schmidt operator on L 2 [0, 2n] ; 

(ii) the eigenfunction corresponding to each non-zero simple eigenvalue of W is a 
2it -periodic solution of (6.3). 

Proof, (i) When a = and A = n 2 with n = 1, 2, . . ., we can take A(x) = sinnx, and 
recover the kernel 

nsinn(x-y) 

W(x,y) = — — — 6.6 

sm(x — y) 

as in the circular ensemble. 

When a^O, there exists by Hochstadt's theorem [18, page 40] an increasing sequence 
(X' n ) which satisfies the estimates 

A^ n _ 1 = (2n-l) 2 + ^ + o(n- 2 ) (n - oo) (6.7) 

0<A / 2n -A' 2n _ 1 =o(n- 2 ), 

and such that (6.3) has a non-trivial solution A. This is Mathieu's function of the first 
kind. As in section 2, we can calculate 

(<9 — *~ ~Q~)^( X '^ = — 2a(sinx cosy + cosx sin y)A(x)A(y), (6.8) 

hence 

W(x, y) = -®J o +V sine a(±(9 + x - V))A(^(6 -x + y))de + g(x - y), (6.9) 

where, by letting x = —y, one finds 

, , A(x/2)A'(-x/2)-A'(x/2)A(-x/2) 

g{x) = : • 

smx 
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Evidently W is a real, symmetrical and continuous kernel, and hence determines a self- 
adjoint and Hilbert-Schmidt operator on L 2 [0, 2n]. 

(ii) By differentiating (6.8) and recalling the definition of W, one can easily deduce 

that 

/ d 2 d 2 \ 

{9-2 - q^j w ( x ' y) = «( cos 2x - cos 2 v) w ( x > y)- ( 6 - 10 ) 



. dx 2 dy' 2 

For v 7^ 0, any non-zero solution / G L 2 [0 7 2tc] of the integral equation 



uf(x)= [ W(x,y)f(y)dy (6.11) 
Jo 

extends to define a twice continuously differentiable and 27r-periodic function on R. Now 
g(x) = f"(x) + acos2xf(x) also gives a 2tt periodic and continuous solution of (6.11); this 
follows from (6.10) by an integration-by-parts argument. By simplicity of the eigenvalue, 
we deduce that g is a constant multiple of /, and hence that / is a 2tt periodic solution of 
Mathieu's equation. 

□ 

Conversely, let Ma be the space of potentials q that have periodic spectrum equal to 
a given A. McKean, van Moerbeke and Trubowitz [19, 20] have shown that Ma can be 
considered as a torus 

M A = {^(A(^) + yj A( Xj ) 2 - 4)°° =i : A 2J _! < x 3 < X 2j ;j = 1,2,...} (6.12) 

over the product over the intervals of instability (A 2 j-i, A 2j ) where A(A) 2 < 4 and that Ma 
is associated with the Jacobi manifold over the Riemann surface of a/A 2 (A) — 4. Hence 
M\ can have dimension n = 0, 1, . . . , oo, equal to the number of simple zeros A(A) 2 — 4. 
The periodic spectrum A is preserved by Hamiltonian flows; in particular, there is a 2ir 
periodic Korteweg-de Vries flow on Ma associated with 

dq dq 1 d 3 q 

d-t =3q dx--2dx^- (03) 

By Theorem 6.1, the potential acos2x gives an infinite-dimensional Ma on which there 
are solutions to KdV that are 2rr periodic in the space variable and almost periodic in 
time [4, Appendix]. This makes an interesting contrast with the concentric KdV equation 
ut + u/(2t) — 6uu x + u xxx = which is used via T 2 x ^ in Theorem 2.3 to linearize Pn for 
the soft edge ensemble [10, page 173]. 
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